Dynamical mean field theory of an effective three-band model for Naa;Co02 
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We derive an effective Ilamiltonian for liiglily correlated t2g states centered at the Co sites of 
Na2,Co02. The essential ingredients of the model are an O mediated hopping, a trigonal crystal- 
field splitting, and on-site effective interactions derived from the exact solution of a multi-orbital 
model in a CoOe cluster, with parameters determined previously. The effective model is solved by 
dynamical mean-field theory (DMFT). We obtain a Fermi surface (FS) and electronic dispersion 
that agrees well with angle-resolved photoemission spectra (ARPES). Our results also elucidate the 
origin of the "sinking-pockets" in different doping regimes. 

PACS numbers: 71.27.-|-a, 71.18.-|-y, 74.25.Jb, 74.70.-b 



The construction of the appropriate low-energy Hamil- 
tonian to describe a highly correlated system is a crucial 
task for an advance in its physical understanding. A clear 
example is the case of the superconducting cuprates. The 
starting point for the description of those materials is a 
three-band model containing the most relevant Cu and 
O orbitals. The parameters of that model were deter- 
mined by constrained-density- functional theory [l[. On 
the basis of the exact solution of the multi-band model 
in a Cu04 cluster (containing one Cu atom and its four 
nearest neighbors) , Zhang and Rice suggested that the es- 
sential low-energy physics of the model is captured by a 
one-band model containing only effective Cu orbitals ^] . 
This has been confirmed by systematic derivations of the 
ensuing one-band Hubbard and t—J models [3|,|3|- These 
models have led to a considerable progress in the un- 
derstanding of the high- Tc cuprates. Similar low-energy 
effective models were derived and used successfully to 
explain the properties of nickelates [5., ^] and other tran- 
sition metal oxides 

In the cobaltates Naa;Co02 a consensus has not yet 
been reached on the appropriate low-energy effective 
Hamiltonian, as different approaches have provided con- 
flicting results. The cobaltates present a clear cut ex- 
ample of strong correlation effects. Not only by its rich 
phase diagram that includes a charge ordered insulator 
and a superconducting state at intermediate dopings, but 
also by the complete failure of the ab-initio band struc- 
ture calculations to describe the shape and topology of 
the Fermi surface measured in ARPES experiments [8|,Q. 
Specifically, first-principles calculations done in the local- 
density approximation (LDA) [lo| predicted a Fermi sur- 
face with six prominent hole pockets along the T — K 
direction, which were never detected in photoemission. 
In addition, the ARPES experiments have revealed the 
presence of dispersive features at the momenta positions 
where those pockets were expected but they were ob- 
served at about 0.2 eV beneath the Fermi surface. Thus, 



they were termed " sinking pockets" and are still awaiting 
a clear physical interpretation. 

Initial theoretical progress was seemingly achieved by 
Zhou et al. who included correlation effects on top 
of a tight-binding model fit to the band-structure from 
first principles calculations in the local-density approxi- 
mation LDA They showed that correlation effects 
may in fact wipe out the pockets by a reduction of the 
bandwidth of the bands crossing the Fermi energy. How- 
ever, the approach of Zhou et al. relied on a simplified 
static Gutzwiller approximation (GA) where the rather 
unrealistic assumption of an infinite strength for the local 
effective t2g Coulomb repulsion U is made. In a differ- 
ent approach to the problem, Ishida et al. [l^ used the 
more elaborate DMFT methodology to treat the correla- 
tion effects, on top of a similar LDA-derived tight-binding 
Hamiltonian. Significantly, the DMFT method allowed 
for the assumption of finite effective Coulomb interac- 
tions. The main finding of that work was the prediction 
that the effect of U is, in marked contrast to the GA, to 
actually increase the size of the LDA pockets, that get 
stabilized due to charge transfer between the and aig 
bands. 

Marianetti et al. 13] using a DMFT calculation simi- 
lar to Ishida et al. found that the pockets can be made 
to disappear for sufficiently large values of U (above 6 
eV), which explained their absence in the infinite U cal- 
culation. Although for realistic values of U the pockets 
still remained, those authors also pointed out that using 
the e'g — aig crystal field splitting as a free fitting param- 
eter, they could eventually be made to disappear. More 
recently, Liebsch and Ishida critically discussed the 
various previous approaches that were based on the LDA 
band-structure as the starting point for the calculation of 
correlation effects. They concluded that, at values C/ ~ 3 
eV which they considered realistic, the presence of pock- 
ets in the Fermi surface is always predicted. They argued 
that this feature, which is in confiict with ARPES data. 
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is robust with respect to the details of the LDA-fits and 
to the form of the interaction term. 

From a more general perspective, one may expect that 
a simple-minded identification of the LDA-derived con- 
duction bands as the relevant manifold where correlations 
are to be included through Hubbard-like interactions in 
an LDA-I-DMFT treatment may not be fully justified 
when systems have a strong covalent character, as it is 
the case of the cobaltates. In particular it was shown 
that the above procedure fails in NiO, and agreement 
with experiments in LDA-I-DMFT calculations is only 
achieved once the O bands are explicitly included in the 
model Interestingly, the results of that approach 

also agree with results from effective models where the O 
atoms have been integrated out using low-energy reduc- 
tion procedures that take into account correlations from 
the beginning 0, [l^ . 

We propose to address the problem of the low-energy 
description of the band-structure of the cobaltates by 
taking a different approach and altogether leave the LDA 
as the starting point of our calculation. Thus, in this 
Letter we perform a low-energy reduction to derive an 
effective Hamiltonian HcS, that includes the determina- 
tion of the values of effective local repulsive interactions, 
and then study its physical behaviour using DMFT. The 
derivation of Hes follows the ideas of previous research 
in the cuprates which used the cell-perturbation method 
[1| and non-orthogonal Zhang-Rice singlets 0, 0] con- 
structions. Basically, the procedure is to divide the sys- 
tem in different cells that are solved exactly, and retain 
their lowest energy states. Then, one includes the inter- 
cell terms along with the effect of the other states as 
perturbations to this low-energy subspace. The result- 
ing effective Hamiltonian differs substantially from those 
previously adopted. In particular, our calculated value 
of U is significantly smaller, again raising questions on 
the justification of the assumption of an infinite value 
for the Coulomb interaction made in Gutzwiller-type of 
approaches. This observation also applies to a recent 
Gutzwiller Density-Functional calculation that reports 
good agreement with ARPES data, in which [/ > 3 — 5 
eV was assumed [l^. A previously derived Hcs also as- 
sumed infinite effective on-site t2g Coulomb repulsions 
instead of calculating them . 

We start from the exact solution of a CoOg cluster 
model containing all 3d orbitals of a Co atom and all 
2p orbitals of its six nearest-neighbor O atoms, assum- 
ing cubic (0;i) symmetry and neglecting spin-orbit cou- 
pling. All interactions inside the 3d shell are included 
[l8|. The parameters were determined fitting XAS ex- 
periments and its polarization dependence [l8|. The re- 
sults, which agree with previous similar studies [l9'|, show 
a large Co-0 covalency and an intra-orbital repulsion 
Urn — 4.5 eV, larger than the Co-0 charge-transfer en- 
ergy. The subscript m refers to the original multiband 
model, to distinguish Um from the corresponding repul- 



sion U of the effective model, which as shown below, is 
strongly reduced due to Co-0 covalency. We recall that 
in the cuprates, Um ~ 10 eV [H, while in their effective 
low-energy one-band Hubbard model [/ ~ 3 eV Q ■ 

The effective model H^s is obtained mapping the 
ground state of the CoOe cluster with four holes onto 
the on-site vacuum of iJoff (no t2g holes at a Co site, 
i.e. Co+'^), and the 6-fold degenerate (spin doublet and 
orbital triplet) ground state for five holes onto the cor- 
responding states with one t2g hole of Hcs- Details of 



the mapping are given in Ref. [17| . We remark that 
Co^'^ and Co+^ in Hcs actually represent highly corre- 
lated states with a Co valence near 2.04 and 2.56 respec- 
tively Hcs reads 



-f^cff — Ho + Hi, (1) 
where the "non-interacting" part can be written as 



Vaa' ^ '-a "a 



(2) 



where dj^^ creates a hole in the t2g orbital a [xy, yz or 
zx) with spin a at site i. However, physically this oper- 
ator represents a non trivial excitation of the same sym- 
metry, which involves also 3d Cg orbitals of Co and 2p or- 
bitals of nearest-neighbor O sites, t' and t correspond to 
the direct Co-Co hopping and to that mediated by O 2p7r 
orbitals [2^ respectively. The latter is the most impor- 
tant one and has been calculated before using many-body 
eigenstates of the CoOg cluster 17[. Finally D accounts 
for the trigonal crystal-field splitting A = 31? between 
e'g and aig orbitals. We take it from quantum-chemistry 
configuration- interaction calculations . These are the 
most reliable methods to determine crystal-field excita- 
tions. Incidentally, it is known that while the LDA may 
provide a good description of the ground state, it does 
not get the energy of excited states right. Therefore it is 
not expected to provide accurate values for D in a highly 
correlated system. Note that although iJg has the form 
of a non-interacting Hamiltonian, the derivation of its pa- 
rameters already involve many-body calculations 17|, |21 1 . 
In fact, similarly as in the studies of cuprates, most of 
the original Co on-site interaction is already included in 
the derivation of HcS through the diagonalization of the 
CoOg cluster. 

The interacting part of i/cft is 
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where the interaction parameters were calculated from 
the comparison between the energy of adding two holes in 
the same CoOg cluster with given symmetry and spin, or 
in different clusters. The eigenvalues of Hj with two holes 
should coincide with the corresponding lowest energy lev- 
els for 6 holes in the cluster. The resulting parameters of 
the model become t = 0.10 eV, D = 0.105 eV, U = 1.86 
eV, U' = 1.27 eV, J = 0.35 eV and J' = 0.17 eV. To our 
knowledge, this the first time that a calculation of the 
interaction terms is reported in this system. Note that 
the values of the U parameters are smaller than those 
used in previous calculations [ll|, 12, [l^l, but are still 
much larger than the bandwidth. We have also added 
a direct hopping between orbitals of the same symmetry 
t' = 0.02 eV that provide the best agreement with exper- 
iments, however our main results are not affected by the 
specific chosen value. 

Here we solve H^s using the DMFT The as- 

sociated quantum impurity problem is a three-orbital 
Anderson impurity that is solved using the Hirsch-Fye 
quantum Monte Carlo (QMC) algorithm [2^. Due to 
the symmetry of the band structure of Hq, the DMFT 
quantum impurity problem and its corresponding self- 
consistency constraint (Dyson Equation) are diagonal in 
orbital and spin indexes. Thus, the resulting local self- 
energies Y^^^^'^ are also diagonal. In order to obtain the 
momentum and energy resolved Green's functions, the 
local self-energies have to be analytically continued to 
the real frequency domain. Thus, we obtained high qual- 
ity QMC data using over one million sweeps to reliably 
perform the continuation by means of a standard maxi- 
mum entropy method |.24j . In the calculations presented 
here, J' and the spin flip terms in Eq. ([3]) were neglected. 
This simplification introduces tiny modifications in the 
results 14]. Thus, we adopt the interaction parameters 
Ua,a = J7 = 1.86 eV for the intra-orbital repulsion, and 
[/^J-f ^U' = 1.27 eV and U^]^, = U' - J = 0.92 eV, 
for the inter-orbital repulsions with opposite or the same 
spin, respectively. 

The predicted band structure is then obtained from 
the imaginary part of the lattice Green's functions given 
by (we study paramagnetic solutions so we drop the spin 
index) Ga(k,u;) = [lo ~ ek,« - 11^^^'^ (lo)]-^ and the 
Fermi surface is mapped out from the lo = crossings of 
the interacting bands. 

We focus our study on the cases that are experimen- 
tally most relevant, i.e., for doping a;=0.3, 0.5 and 0.7, 
which range from stronger to weaker correlations. For 
reasons of space, the data displayed in the figures are for 
a;=0.3 and 0.7. Due to their high computational cost, the 
lowest temperature that we study is T = 360 K. Compar- 
ison with calculations at higher temperatures 720 K) 
indicates that we have indeed achieved the low T limit. In 
addition, as will be shown latter, the width of the quasi- 
particle band at the Fermi energy, which is the smallest 
energy scale in the electronic structure, is much larger 



than the temperature of the calculation. 

In Fig. [1] we show our results for the evolution of the 
Fermi surface as function of increasing doping along with 
the respective experimental ARPES data. We observe 
good agreement in the shape and size of the FS at all 
doping levels. Significantly, the hole pockets are absent 
in our results. 

The experimental FS for x = 0.3 is somewhat more 
rounded than the theoretical one. This may partialy be 
due to the relatively large thermal broadening in the cal- 
culation, but may also be due to lack of hopping terms 
at longer distances, beyond those included in Hq. 




X = 0.72 
i= 0.60 
X = 0.48 
0.0 • J! = 0.30 
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FIG. 1: (Color online) (a) Calculated FS for dopings x = 0.3 
and X — 0.7 (inner (hole) edge of larger and smaller hexagon 
respectively), (b) ARPES measurements from 

The details of the band structure are shown in Fig. [2l 
We observe that the data reveal several contributions 
that can be associated to either coherent (i.e. quasipar- 
ticle like) or incoherent (i.e. Hubbard) bands. The in- 
coherent bands are characterized by dispersive structures 
similar to those of the "non-interacting" Hamiltonian Hq 
(though usually less defined due to shorter lifetimes) that 
appear far from the Fermi energy. These large energy 
shifts are of course due to the local interactions of Hes 
[Eq. Q]. In the top left panel of Fig. [2] we show the 
full band-structure for the strongly correlated case x = 
0.3. There, one can observe several incoherent bands 
that appear shifted down in energy, at ^ -1, -1.75 and 
-3 eV, with the first one earring a large part of the spec- 
tral intensity. Their shapes reveal their dominant orbital 
content and their energy shifts can be understood from 
the values of the inter- and intra-orbital Coulomb repul- 
sions. At higher dopings, the correlation eflFects decrease 
and these incoherent bands rapidly loose spectral inten- 
sity. On the other hand, the coherent bands are near 
the Fermi energy and their band structure is somewhat 
narrowed with respect to that of Hq [Eq. ([2])] due to the 
effect of Hi [Eq. ([3])], indicating the enhancement of the 
effective mass. 

The top right panel of the figure shows details of the 
band structure at a;=0.3 and, for comparison with the less 
correlated case, the lower right panel shows similar data 
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for x=0.7. Interestingly, these results reveal a novel in- 
sight on the nature of the "sinking pockets" , whose exper- 
imental data we reproduced in the lower left panel. We 
find that while the sinking pockets are present at both, 
low and high dopings (they are indicated by boxes in the 
respective panels), their physical origin is qualitatively 
different. At higher x, correlations are low and the band 
structure does not differ much from the non-interacting 
case. Thus, the sinking pocket in this case can be simply 
associated to the top of the band with mostly e'g char- 
acter [indicated by a box in panel (d)]. In contrast, at 
X — 0.3, in the strong correlation case, as we discussed 
before, the band structure is dramatically modified and 
that interpretation is no longer possible. In fact, the 
strongest contribution to the e'g band is shifted down in 
energy by about 1 eV. This shift is due to the inter-orbital 
Hubbard repulsion, and can be more easily understood 
in a hole picture. As there is about one hole in the aig 
band, putting a second hole costs Uaig,e'g ^ 0.92 — 1.27 
eV if the hole goes into the e'g band [or Uaig,aig ^ 1.86 
eV into the aig band [see panel (a)]. However, in the 
ground-state there is also a non-negligible amplitude for 
a configuration with no holes in the aig band, thus one 
may create a hole in the original e'g band with no extra 
Coulomb energy cost. Such a state would have a re- 
duced spectral intensity but a similar dispersion as that 
of the non-interacting e'g band, thus accounting for the 
sinking pocket. This is confirmed by our calculations on 
the weight of the states with different symmetries in the 
coherent bands (not shown). A clear signature of this 
sinking pocket state is indicated by a box in our numer- 
ical data of panel (b). We note that our results show 
a lower energy edge at ~ —0.2 eV in good agreement 
with aU available ARPES data [1, Nevertheless, the 
experimental situation is less clear for the determination 
of the dispersive shape of the sinking pockets at higher 
binding energies. With regard to the influence of possi- 
ble Na ordering, recent ARPES experiments |9|] conclude 
that that feature does not affect significantly the position 
of the sinking pockets in the unfolded Brillouin zone. 

In conclusion, motivated by the apparent failure of 
LDA band structure to provide a sound starting point for 
the calculation of strong correlation effects in the cobal- 
tates, we derive an effective low-energy Hamiltonian that 
includes the strength of the local Coulomb repulsive in- 
teractions. The effective model is obtained from finite 
cluster and quantum chemistry calculations, with essen- 
tially no adjustable parameters. The effective Hamilto- 
nian is treated with DMFT to compute the effects of 
correlations at different doping levels. We find that the 
evolution of the Fermi surface is in good agreement with 
the experimental data. Importantly, the LDA-predicted 
hole pockets, that are not seen in the ARPES data, are 
also not present in our results. One difference with re- 
spect to LDA calculations is how the effect of correlations 
separates the bands as discussed above. However, for re- 




FIG. 2: (Color online) (a) Full band structure for x = 0.3. 
The X = 0.3 and x = 0.7 cases are detailed in (b) and (d) re- 
spectively, where the structures in the dispersion correspond- 
ing to the sinking pockets are highlighted by a box. (c) sinking 
pockets for various dopings, with comparison to LDA [^. 



alistic U this feature alone is not enough to destroy the 
pockets and the value of D, usually underestimated by 
LDA, was also shown to play a key role. We obtained 
the detailed interacting electronic structure that reveals 
sinking pockets at all dopings. Significantly, their origin 
is qualitatively different in the high and low doping cases. 
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